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Resume : La cartographic de I'univers en radio, a travers I'observation de la raie 
a 21 cm de I'hydrogene atomique, constitue une approche complementaire aux re- 
leves optiques pour 1' etude des grandes structures, et en particulier des oscillations 
acoustiques baryoniques (BAO). Nous proposons une methode originale de la me- 
sure de la distribution de I'hydrogene atomique neutre a travers une cartographic a 3 
dimensions de remission du gaz a 21 cm, sans passer par 1' identification des sources 
compactes (galaxies . . . )• Cette methode necessite un instrument a grande sensibi- 
lite et ay ant une grande largeur de bande instantanee (> 100 MHz), alors qu'elle 
pent se contenter d'une resolution angulaire moyenne (10 arcmin). L'instrument de- 
vra avoir une surface de collection de quelques milliers de m^ et quelques centaines 
de lobes simultanes. Ces contraintes peuvent etre satisfaites avec un reseau dense 
de recepteurs en mode interferometrique ou un reseau phase au foyer d'une grande 
antenne. 

1 Introduction 

Les questions concernant la nature et les proprietes de I'energie noire se trouvent au coeur de la 
cosmologie et de la physique aujourd'hui. L' energie noire est une composante mysterieuse responsable 
de I'acceleration de I'expansion de I'Univers. Les oscillations acoustiques du plasma photons-baryons 
anterieur au decouplage sont a I'origine des modulations du spectre de puissance de la distribution de 
matiere, appelees BAO. Ces oscillations acoustiques ont clairement ete mises en evidence dans les ani- 
sotropics du fond diffus Q. La mesure de I'echelle de longueur des oscillations a differents redshifts est 
une sonde cosmologique de type regie standard, consideree comme I'une des methodes les plus robustes 
pour contraindre I'equation d'etat de I'energie noire [2J. 

La distribution de matiere dans I'univers peut bien sur etre observee en optique, en utilisant les ga- 
laxies comme traceur. L' observation radio de remission a 21 cm de I'hydrogene neutre est une methode 
complementaire pour determiner la distribution de matiere. U-L Pen et J. B. Peterson ont propose un 
projet radio en vue de cartographier la distribution de I'hydrogene atomique a travers son emission a 21 
cm (transition hyperfine, 1.42 GHz a z=0). II est important de noter que la methode d' observation envi- 
sagee ici repose sur la cartographic de remission totale a 21 cm en fonction du decalage vers le rouge, 
sans rechercher I'observation des sources individuelles (galaxies. . .) iSl O. L'objectif est de produire 
une cartographic a 3 dimensions de remission HI (2 angles + frequence ( redshift)), avec une resolution 
angulaire d'une dizaine de minutes d'arc et de quelques dizaines de kHz en frequence (dz/z < 10~^). 

Un programme de recherche et developpement en electronique pour I'interferometrie radio a ete 
initie par le LAL (IN2P3) et I'Irfu (CEA) fin 2006, auquel participent egalement 1' observatoire de Paris 
(GEPI et station de Nangay) depuis I'ete 2007. Une chaine complete comprenant des modules analo- 
giques de conditionnement des signaux radio, des modules de numerisation et de filtrage numerique et 
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un sy Sterne d' acquisition et de traitement haute performance, a ete developpee grace au soutien de P2I, 
de rirfu (CEA), de l'IN2P3 (CNRS) et du PNCG. Cette chaine est utilisee pour caracteriser les proto- 
types de reflecteurs cylindriques developpes par I'universite Carnegie-Mellon (CMU) a Pittsburgh. Le 
projet FAN (Focal Array at Nangay) utilise egalement la chaine electronique et informatique developpee 
pour BAORadio. Ce projet a pour but le developpement et la caracterisation d'un prototype de reseau 
de recepteurs au foyer du NRT (reseau phase) et pourrait deboucher a terme sur un systeme multi lobes 
qui augmenterait de plus d'un ordre de grandeur la sensibilite du radio telescope de Nangay en mode 
cartographic iSJ. 

Les enjeux scientifiques et des difficultes specifiques de la methode d' observation envisagee, la carto- 
graphic 3D a travers la mesure de la temperature de brillance en fonction de la frequence (Total Intensity 
Mapping) seront presentes dans le paragraphe[2j Nous discuterons en particulier le niveau de sensibilite 
necessaire pour la mesure du spectre des inhomogeneites de densite P(k) et les difficultes de separation 
du signal cosmologique des avant-plans radio. Les developpements techniques qui ont ete menes de- 
puis 2007 seront presentes dans le paragraphe [3j Les resultats obtenus lors des tests de qualification 
de r electronique a Nangay (NRT) et sur les prototypes de reflecteurs cylindriques de Pittsburgh seront 
brievement discutes. 



2 Cartographic 3D de remission a 21 cm 

L' etude des proprietes statistiques de la distribution de matiere dans I'univers permet de tester la 
validite du modele cosmologique et pent etre utilise pour determiner les valeurs des parametres cosmo- 
logiques, en particulier celles des densites cosmiques. On determine la distribution cosmique de matiere, 
dominee par la matiere noire en utilisant des traceurs, les galaxies dans la majorite des releves. Ces ga- 
laxies sont en general detectees a travers leur emission optique, la position etant alors determinee par 
imagerie optique, et le decalage vers le rouge {z) par spectroscopic. Une methode similaire pent etre uti- 
lisee dans le cas de I'observation des galaxies et des nuages compact de gaz en radio, a travers remission 
a 21 cm. Le decalage vers le rouge de la source est alors determine directement en comparant la frequence 
de I'onde radio observee v a celle de la raie a 21 cm {z = 1420MHz _ applications cosmologiques 

de I'observation des galaxies a 21 cm sont par exemple discutees dans la reference (61. 

Mais la faible luminosite des galaxies a 21 cm limite la possibilite de detection de ces objets au 
voisinage de notre Galaxie avec les instruments actuels. Une antenne d'une surface de A ^ lO^m^ 
comme celle du grand radiotelescope de Nangay (NRT) ne permet la detection des galaxies a 21 cm 
{Mhj ^ lO^^M©) que jusqu'a z < 0.1 avec un temps d'integration d'une heure {Sum ^ ImJy). 
La detection des galaxies en radio a 21 cm a des distances cosmologiques (z > 1) ne peut done etre 
d'envisagee qu'avec des surfaces de I'ordre de A ^ lO^m^. 

Le spectre des inhomogeneites de densite P{k) est I'un des outils mathematiques les plus utilises 
pour quantifier et analyser les proprietes de la distribution de matiere. Si on note p{f) la densite de 
matiere dans le cosmos et F(fc) la transforme de Fourier du champ des inhomogeneites de densite S{f) : 

5(f) = ^ - 1 ^TF^ F{k) 
P 

Si on admet I'isotropie a grande echelle de la distribution de matiere, les proprietes statistiques de la 
distribution de matiere peuvent etre resumees dans le spectre de puissance P{k) 

P{k) =< \F{k)\^ > avec k = \k\ 

Les Oscillations Acoustiques Baryoniques (BAO) correspondent aux empreintes laissees dans la dis- 
tribution de matiere par les ondes de pression (son) qui se propageaient dans le plasma photon-baryon 
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PNoise(k) : 5/10/100 beams/polar @z=1 
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Figure 1 - Sensibilite d'un instrument radio pour la mesure du spectre des fluctuations de temperature 
de briflance de remission Hi a z = 1. Observation du quart du ciel (tt steradians) pendant une annee. 
Trois configurations instrumentales, avec un niveau de bruit Tgys = 50 K : (a) antenne de D=75 metres 
de diametre, equipee d'un instrument focal (FPA) avec 5 lobes, (b) D=100 metres et 10 lobes, (c) D=100 
metres et 100 lobes Q. 

avant le decouplage. Les BAO se manifestent conmie une legere modulation (< 10%) du spectre de 
puissance P{k) pour kn ^ isolXpc • Cette echelle de longueur (150 Mpc) correspond des structures 
dans le plan transverse ayant une taille angulaire de ^ 1 degre. La temperature moyenne de brillance 
a 21 cm et le niveau des fluctuations attendues sont inferieurs au milli-Kelvin. Un instrument radio 
de taille modeste (D < 100m) permettrait done de resoudre ces structures. Par contre, il faut un ins- 
trument ayant une grande sensibilite a la temperature de brillance pour eviter que le signal cosmolo- 
gique soit noye par les fluctuations du bruit instrumental, quantifie par la temperature systeme (Tgys)- 
La determination du spectre de puissance necessite en effet le sondage d'un grand volume de I'uni- 
vers afin de limiter la dispersion due au nombre de modes mesure (Sample variance) . II faut done 
realiser un releve couvrant une fraction notable du ciel 10000 deg^) et une tranche en decalage 
vers le rouge Az > 0.2. Ces contraintes peuvent etre satisfaites par un instrument ayant un grand 
champ de vue instantane (10-100 deg^) et a large bande (> 100 — 200MHz). La cartographic a 3 di- 
mensions (2 angles, frequence/redshift) de la temperature de brillance avec une resolution angulaire 
Sft ^ lOarcmin x lOarcmin, et une resolution en frequence Su ^ lOOkHz pourrait etre utilisee pour la 
determination du spectre de puissance de la distribution de I'hydrogene atomique. La figure [TJmontre la 
contribution du bruit instrumental, pour plusieurs configurations instrumentales, avec 5,10 et 100 lobes 
simultanes, comparee au spectre de puissance attendu de la temperature de brillance de remission a 21 
cm a z = 1. 

Le niveau du bruit instrumental ne represente pas le seul defi experimental de la cartographic 3D de 
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Figure 2 - Gauche : temperature de brillance de remission synchrotron et des sources radio continuum 
a 840 MHz (echelle des couleurs en Kelvin), pour une zone de 30 x 30 deg^, parmi les plus froides de 
la Galaxie. Droite : temperature de brillance attendue pour remission a 21 cm de I'hydrogene atomique 
a z=0.7 (echelle des couleurs en milli-kelvin) ^ 



remission a 21 cm. En effet, les avant-plans radio dus a remission synchrotron de la Voie Lactee et des 
radio sources sont plus d'un millier de fois plus intenses que remission de I'hydrogene atomique neutre. 
La figure [2] montre le niveau des variations de temperature de brillance attendu pour la raie a 21 cm {u 2^ 
840MHz, z 0.7) comparees a celles dues au synchrotron galactique et aux sources radio. Neanmoins, 
des etudes preliminaires montrent qu'il devrait etre possible de separer les deux contributions en utilisant 
la forme du spectre en frequence des emetteurs radio Q. En effet, le synchrotron galactique et les radio 
sources ont des spectres d' emission en loi de puissance oc u^. Cette caracteristique pent etre mise a profit 
pour isoler le signal Hi cosmologique. 

3 Uelectronique BAORadio 

L'originalite du systeme electronique developpe reside dans sa conception entierement numerique, et 
dans r utilisation de circuits cables a base de FPGA afin d'obtenir une puissance de traitement importante, 
avec un cout et une consommation moderes. La gamme de frequence 0.5-1.5 GHz est decoupee en plu- 
sieurs bandes de 250 MHz de large, et le signal analogique provenant des recepteurs est numerise a 500 
MHz, apres une etape d' amplification et de filtrage. Ce flot numerique est ensuite traite par les differents 
etages de la chaine electronique et d' acquisition. D'autres developpements de systemes numeriques pour 
des instruments radio |9 1 ont ete realises ou sont en cours : CASPER/ROACH a Berkeley LIOJ et Uni- 
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board en Europe |TT|. 

Des elements de la chaine electronique et le systeme d' acquisition ont ete congus et realises durant 
les trois dernieres annees (2007 a 2009). Outre des tests de mise au point en laboratoire, le programme 
de travail comprend des tests de qualification effectues en situation reelle aupres du radiotelescope de 
Nangay et du prototype de Pittsburgh. 

1. Module analogique d' amplification, de filtrage et de decalage en frequence. Ces modules, congus 
et realises par I'lRFU, permettent le decoupage du signal en bande de 250 MHz. 

2. Systeme de distribution d'horloge et de signaux de synchronisation (DCLK). Celui-ci, egalement 
sous la responsabilite de I'lRFU, permet d'assurer le synchronisme temporel avec la precision 
necessaire. 

3. Une carte de numerisation a 4 voies, capable de numeriser les signaux a 500 MHz a ete congue 
et realisee par le LAL. Cette carte est equipee de circuits programmables puissants, capable d'ef- 
fectuer au vol le filtrage numerique des signaux pour les separer en composantes quasi monochro- 
matiques (FFT). Les donnees brutes ou apres filtrage (FFT) sont transmises par des liens optiques 
haut debit (5 Gbits/s) vers les ordinateurs d' acquisition, ou vers un systeme dedie de correlateur 
et de reconstruction de lobe. Le code embarque (firmware) de filtrage numerique est realise par 
riRFU et le LAL. 

4. Un module de reception des donnees, a la norme PCI-Express, capable de soutenir des debits de 
plusieurs centaines de mega-octets par seconde permet le transfert des donnees vers la memoire 
des ordinateurs d' acquisition. Le code embarque de ce module est egalement developpe au LAL. 

5. Un systeme d' acquisition et de traitement au vol des donnees sur des grappes de PC a ete developpe 
sous la responsabilite du LAL, en collaboration avec I'lRFU. 

6. La conception et la realisation d'un prototype d'un systeme de correlateur/synthetiseur de lobes 
d'antenne a base de FPGA est en cours a I'lRFU. Les premiers tests de ce prototype sur le ciel 
devront avoir lieu dans les prochains mois. 

Les premiers tests de la chaine electronique ont ete effectues avec succes a Nangay au foyer du ra- 
diotelescope decimetrique (NRT) en juillet 2008. D'autres campagnes de tests ont ete menees durant 
I'hiver et le printemps 2009. lis ont permis en particulier de tester des versions ameliorees des codes 
embarques (firmware FFT) et de 1' acquisition multi-voies. Plusieurs campagnes de mesure en mode 
interferometrique ont ete effectuees sur le prototype des reflecteurs cylindriques de Pittsburgh en Juin 
2009, novembre 2009 et decembre 2010. Un programme d' observation des amas de galaxies utilisant 
les recepteurs cryogenique du NRT et I'electronique BAORadio, ainsi que les mesures avec le prototype 
FAN ont debute fin 2010 a Nangay. 

Les observations en mode interferometrique realisees sur le prototype de PittsBurgh en Novembre 
2009 ont permis de verifier le fonctionnement de la chaine electronique : 1' observation de sources 
brillantes (CasA, CygA, Soleil) a permis la calibration des voies d' electronique, la reconstitution de 
lobes numeriques (de type reseau phase) ainsi que la verification de la stabilite temporelle relative des 
voies d' electronique. Une partie des resultats obtenus est resumee sur la figure |4] 

4 Conclusions, perspectives 

La cartographic a 3 dimensions de la brillance en Hj par la methode de la mesure de la temperature 
de brillance est une methode prometteuse pour des releves cosmologiques de la distribution de matiere 
jusqu'a z < 3, et pour contraindre I'energie noire en utilisant les Oscillations Acoustiques Baryoniques 
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Figure 3 - Carte de numerisation (4 voies a 500 MHz) et de filtrage numerique des signaux developpee 
pour le projet BAORadio. 
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Figure 4 - Signaux d' interference et synthese de lobe, donnees Pittsburgh de Novembre 2009 



6 



(BAO). Un instrument multilobes (FPA ou interferometre) ayant un grand champ de vue instantane (10- 
100 deg^) et une bande large en frequence (quelques centaines de MHz) devra etre utilise afin d'obtenir 
une sensibilite suffisante pour le sondage d'un grand volume de I'univers a 21 cm. La chaine electronique 
presentee ici a ete congue pour equiper un systeme multi-recepteurs, comprenant plusieurs centaines de 
voies de numerisation. L' electronique analogique et les cartes de numerisation peuvent etre reparties 
sur une surface de plusieurs hectares, les signaux numeriques etant achemines par fibre optique vers un 
systeme central de traitement de donnees au vol. Ce systeme de traitement comprendra une grappe de 
calculateurs relies par un reseau haut-debit, equipes de cartes de traitement a base de GPU ou FPGA. Le 
cout de la chaine electronique et du systeme central de traitement de donnees devrait se situer dans une 
fourchette de 3000 a 5000 euros par voie de numerisation. 
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